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CHAPTER  1 


INTRODUCTION 


The  physics  of  laser  induced  sound  in  liquids  is  divided  into  two  parts, 
the  coupling  of  the  electromagnetic  energy  of  the  light  field  into  the  fluid  to 
create  a  disturbed  region,  from  which  acoustic  energy  is  radiated,  and  the  pro¬ 
pagation,  subsequently,  of  the  radiated  acoustic  wave.  In  this  report,  I  present 
a  brief  examination  of  a  nonlinear  aspect  of  the  first  of  these  problems,  for  the 
case  of  surface  absorption. 

Broadly,  two  principal  coupling  mechanisms  of  laser  energy  are  absorption 
by  particulate  impurities  in  the  liquid,  and  absorption  by  dielectric  loss. 
Depending  upon  laser  intensities  and  wavelength,  dielectric  loss  in  the  fluid 
may  also  entail  dielectric  breakdown.  Particulate  absorption*  and  breakdown^ 
each  has  been  observed  in  focal  regions  beneath  the  liquid  surface,  for  water 
using  ruby  lasers  (X  =  0.6943  urn).  Miniature  volume  explosions  result  in  the 
production  of  acoustic  energy.  Early  reported  observations  of  laser-induced 
sound,  again  using  a  ruby  laser,  appeared  consistent  with  simple  linear  dielectric 
absorption  below  threshold  for  fluid  breakdown,  leading  to  thermal  generation  of 
an  acoustic  stress  wave.3 

Early  investigations  also  proved  the  utility  of  Nd:glass  (X  =  1.06  um)  and  of 
CO2  (X  =  10.6  pm)  lasers;  in  addition,  dyes  have  been  used  to  control  the  optical 
absorption  constant.4  The  linear  optical  absorption  length, 


a-1  =  i<ni  , 


ni  =  Im  n 


(1.1) 


where  n  is  the  complex  refractive  index  of  the  fluid,  varies  with  wavelength  X 
over  as  many  as  ten  orders  of  magnitude  in  water.  Thus  in  the  visible  region, 
ct~l  varies  in  water  from  a  meter,  or  so,  at  the  red  end  to  several  tens  of  meters 
in  theQblue-green  window,  while  in  the  ultra-violet,  around  600^,  it  falls  as  low 
as  100A(=0.01pm) .  On  the  long  wavelength  side,  for  Nd:glass  laser  light  at 
1.06  um,  ot“l  =  6.0  cm,  and  it  falls  precipitously  again  (though  non-monotonical ly) , 
four  orders  of  magnitude  to  11.8  um,  at  X  =  10.6  um  (CO2) . 


Bell,  C.  E.,  and  Landt,  J.  A.,  Appl .  Phys .  Lett .  ,  Vol .  10,  1967,  p.  46. 

"'Barnes,  P.  A.,  Studies  of  Laser-Induced  Breakdown  Phenomena  in  Water, 

Ph.D.  Thesis,  Simon  Fraser  University,  1969. 

^Carome,  E.  F.,  Moeller,  C.  E.,  and  Clark,  N.  A.,  Appl .  Phys .  Lett . ,  Vol.  4, 
1964,  p.  95.  *  . 


+Gournay,  Luke  S.,  J .  Appl .  Phys . ,  Vol.  40,  1966,  p.  1322. 
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Figure  1  shows  the  behavior  of  the  real  and  imaginary  parts  of  n  as  functions  of 
frequency,  f  =  c/\,  for  water. 5  Note  the  strong  absorption  band  centered  at 
'  =  2.950  cm,  where  a" 1=0.840  pm;  see  Figure  2  taken  from  recent  reported  data 
in  the  infrared. h  From  tabulated  data  in  Ref.  6,  the  refractive  index  of  water 
at  \  =  2.950  pm  is  n  =  1 . 317+0. 2  82 i ,  and  at  \  =  10.6  pm,  n  -  1.18+0.075i. 

For  small  values  of  penetration  of  light  into  the  medium  is  effectively 

blocked,  and  the  region  of  disturbed  fluid  is  confined  to  a  thin  surface  layer, 
in  the  COj  laser  case,  where  i  10“^  cm,  3.10  j/cm^  is  the  minimum  value  for 

t ime- integrated  impulse  laser  intensity  necessary  to  vaporize  one  absorption 
depth  of  water  from  room  temperature  (20°C) .  The  first  reported  CO2  laser- induced 
sound  study  in  water  is  due  to  Bunkin,  et  al^,  who  observed  surface  vaporization. 
(.Vaporization  also  had  been  observed  previously  under  other  conditions  by 
Gournay.^)  Numerical  calculations  were  performed  subsequently  by  Feiock  and 
Goodwin,®  who  used  a  one-diinensional  hydrocode  for  an  equilibrium  vaporization 
model,  which  included  an  inhomogeneous  Beer's  Law  that  assumed  proportionality 
of  a  to  where  0  is  the  ( 1  iquid- vapor)  medium  density. 

In  fact,  surface  vaporization,  for  modest  conditions  of  laser  intensity  and 
time-history,  is  readily  achieved  and  is  explosively  sudden.  A  strong  air  shock, 9 
F/Pq  =  47  bar,  was  formed  above  the  water  surface  by  absorption  of  a  1.67j  CO2  TEA 
laser  pulse  (150  ns  FIVUM)  having  a  film  burn  diameter  0.80  cm.  The  spot  diameter 
on  the  water  surface,  as  defined  by  the  extent  of  maximal  fluid  disturbance 
(shadowgraph  pictures)  was  0.68  cm.  The  corresponding  laser  fluence  (intensity 
impulse)  is  1.67/--  x  0.682  =  4.6j/cm2,  an  upper  limit  estimate  only  slightly 
above  the  threshold  scale  value  of  3.1  j/cm-  cited  earlier  for  vaporization  of 
one  absorption  depth. 

The  weak  coupling  case  of  low-level  laser  input,  where  only  thermal 
expansion  will  occur,  but  no  surface  vaporization,  gives  acoustic  pressures 
proportional  to  K/Cp,  where  K  is  the  volume  coefficient  of  expansion  and  Cp  the 
specific  heat  of  the  fluid.  For  water,  K  =  2 . Ixl0_4 (K° ) " * ,  and  Cp  =  4.2j/g-K°, 
giving  K/Cp  =  4. 7x10" 12(cm/sec) “2 .  For  benzene,  the  same  quantity  has  the  value 
62. 3x10“ ^“tem/sec) which  is  an  order  of  magnitude  larger.  The  physics  of 
laser- induced  thermoacoustic  effects  has  received  extensive  treatment  in  recent 
(Soviet)  literature,  and  has  been  reviewed  recently  by  Lyamshev  and  Sedov. ^ 


Gournav,  Luke  S.,  J .  App 1 .  Phys . ,  Vol.  40,  1966,  p.  1322. 

J. Jackson,  J.  D. ,  Classical  Electrodynamics  (2nd  Ed.)  (New  York:  John  Wiley  & 
Sons,  Inc.,  1975),  p.  291. 

^Downing,  Harry  D. ,  and  Williams,  Dudley,  J.  Geophys.  Res.,  Vol.  80,  1975, 
p.  1956. 

^Bunkin,  F.  V. ,  Karlov,  N.  V.,  and  Komissarov,  V.  M. ,  Sov.  Phys. — JETP  Lett., 
Vol .  13,  1971,  p.  341. 

'"Veiock,  F.  D.,  and  Goodwin,  L.  K.  ,  J .  Appl .  Phys ,  ,  Vol.  43,  1972,  p.  5061. 
^Be 1 1 ,  C.  E. ,  and  Maccabee,  B.  S.,  Applied  Optics,  Vol.  13,  1974,  p.  605. 
^Lyamshev,  L.  M. ,  and  Sedov,  L.  V.,  Sov.  Phys. — Acoustics,  Vol  27,  1981, 


FREQUENCY  (Hz» 

NOTE:  THE  VISIBLE  REGION  OF  THE  FREQUFNCY  SPECTRUM  IS  INDICATED  BY  THE  VERTICAL 
DASHED  LINES.  THE  ABSORPTION  COEFFICIENT  FOR  SEA  WATER  IS  INDICATED  BY  THE 
DASHED  DIAGONAL  LINE  AT  THE  LEFT. 

FIGURE  1.  THE  ABSORPTION  COEFFICIENT  FOR  LIQUID  WATER  AS  A  FUNCTION  OF  LINEAR  FREQUENCY 
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in  contrast,  theories  of  the  surface  vaporization  processes,  also  reviewed 
very  recently,  by  Lyamshev  and  Naugol*  nykh'11  --  who  include  the  topics  of 
nonlinear  thermoacoustics  and  optically  induced  dielectric  breakdown,  as 
well  --  are  much  less  well-developed.  Two  cases  (at  least)  have  to  be  distin¬ 
guished:  (1)  High  laser  intensity,  where  a  plasma  is  formed  above  the  surface. 

In  the  model  of  Pirri,^  sensibly  all  of  the  energy  absorption  occurs  in  the 
plasma,  which  forms  the  leading  edge  of  a  detonation  front  (laser  supported 
detonation  wave)  that  advances  in  the  opposite  direction  to  that  of  the  laser 
beam  propagation.  The  force  delivered  to  the  water  surface  is  the  result  of 
shocking  the  air  to  high  pressure.  (2)  Low  laser  intensity,  where  rapid  and 
sustained  vaporization  occurs,  but  no  plasma  front  is  created.  In  this  case,  the 
laser  energy  deposition  goes  directly  into  the  vapor  and  the  liquid.  In  both 
cases,  a  flow  field  must  be  determined  in  order  to  predict  the  full  fluid 
disturbance  from  which  the  acoustic  pulse  is  radiated;^  both  the  absorption  and 
stress  coupling  processes  in  the  two  cases  are  different.  I  will  not  consider 
case  1  here,  but  only  case  2. 

It  is  not  clear  that  the  flow  in  the  vaporization  process  is  steady  since 
the  vaporization  layer  itself  will  be  no  more  than  a  few  molecular  diameters 
thick,  and  is  probably  unstable. 14  The  absorption  properties  of  the  vapor  and 
heated  water  to  the  incident  light  will  vary  with  time  and  with  space  through  the 
density  and  temperature  dependence  of  the  electromagnetic  response  functions 
(dielectric  constant),  even  without  a  surface  instability  of  the  vaporization 
layer.  So,  the  complete  laser  field  will  be  affected  in  the  fluid.  Except  for 
approximate  allowances  for  dielectric  inhomogeneity  (e.g.,  Ref.  ^) ,  this  aspect 
of  the  laser-liquid  coupling,  problem  appears  not  to  have  received  any  attention, 
so  far;  that  aspect  is  the  subject  of  the  present  report.  The  main  purpose  of 
an  investigation  of  this  kind  is  to  search  out  and  identify  mechanisms  for 
optimizing  laser  energy  coupling.  Accordingly,  I  will  not  restrict  the  analysis 
to  the  case  of  10.6  pm  radiation  in  water,  exclusively,  but  will  be  concerned 
more  generally  with  nonlinear  time -dependent  dielectric  coupling  dynamics, 
broadly  appropriate  to  the  surface  evaporation  process. 

The  plan  of  the  report  is  as  follows.  In  Section  2,  I  specify  the  physical 
basis  of  the  equations  to  be  derived  subsequently,  which  is  the  constitutive 
relation  for  adiabatic  non- stationary  response  of  an  inhomogeneous,  local 
dielectric.  In  Section  3,  I  derive  the  source-free  Maxwell  equations,  for  a 
medium  defined  in  the  previous  section.  In  Section  4,  I  develop  the  slow  time- 
scale  approximation  scheme  and  derive  equations  of  the  evelope  fields  of  a  pulse 


« 


a 


^Lyamshev,  L.  M. ,  and  Naugol'nykh,  K,  A.,  Sov.  Phys. — Acoustics,  Vol.  27, 
1981,  p.  337 . 

^Pirri,  Anthony  N.,  Phys.  of  Fluids,  Vol.  16,  1973,  p.  1435. 

^  Lighthill,  M.  J. ,  Proc.  Roy.  Soc .  London,  Vol  A211,  1952,  p.  564,  and 
Vol.  A222,  1954,  p.  1.  See  also  Robert  T.  Beyer,  Nonlinear  Acoustics, 
written  for  the  Naval  Sea  Systems  Command,  1974. 

14  ' 

Anisimov,  S.  I.,  Tribel'skii,  M.  I.,  and  Epel'baum,  Y.  G.,  Sov.  Phys. — 

■  JETP,  Vol.  51,  1930,  p.  802.  ~  . " 
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wave  with  a  high  frequency  carrier.  In  Section  5,  I  specialize  these  equations 
to  the  case  of  a  one-dimensional iy  stratified  medium,  and  appropriate  to  a  wave 
incident  normally  in  the  fluid.  I  also  give  a  short  physical  discussion  of 
limiting  cases,  leaving  the  more  detailed  analyses  needed, for  later.  In  Section 
b,  I  summarize  the  report,  and  make  a  couple  of  concluding  remarks. 


CHAPTER  2 


CONSTITUTIVE  RELATION 


The  complex  dielectric  function  for  water  is  a  function  of  density  p,  and  of 
temperature  T  also.  This  can  be  significant  in  the  surface  regions  especially, 
and  at  the  vaporization  front.  Thermal  and  hydrodynamic  response  to  the  laser 
field  create  spatial  inhomogeneity  and  time-dependence  (non-stationarity)  of 
medium  response  functions  so  that  the  linear  (and  local)*  dielectric  displacement 
is  given  by 


D(x,t)  =  J*dt '  e  (x,  t ;  t-t '  )E  (x,  t ' )  . 


(2.1) 


The  time-scale  for  the  explicit  t-dependence  shown 
of  p  and  T,  and  will  be  assumed  very  long  and  slow 
(carrier)  period  f“l  =  A/c  =  0.03  ps,  for  A  =  10.6 
may  be  regarded  as  adiabatic,  we  may  write 


for  e  is  that  for  variation 
in  comparison  with  the  laser 
am.  If  the  t-dependence  of 


£(x,t;  t-t')  %  e(p,T;  t-t') 


(2.2) 


where  the  quantity  on  the  right  hand  side  is  the  response  function  for  a 
stationary  medium,  having  density  and  temperature  given  by  the  values  in  the 
fluid  at  time  t.  The  explicit  (x,t) -dependence  required  is  then  assumed  to 
enter  through  that  realized  in  the  dynamics  of  the  fluid  for  p  and  T;  viz. 
p=p(t, t)  and  T=T(2,t) .  Eq.  (2.  1)  now  becomes 

D(x,t)  =/dt  'e  (p(x,t)  ,T(x,t)  ;  t-t '  )E(x,t ' )  (2.3a) 


=  eop(o,T)E,  (2.3b) 

where  e0p  denotes  the  integral  operator  in  the  first  equation.  The  limitation  to 
the  validity  of  the  quasi-static  approximation,  p>qs.  (2.2)  and  (2.3),  cannot  be 
assessed  without  a  microscopic  (molecular  radiator)  theory  of  the  laser  field 
absorption  process,  which  is  outside  the  scope  of  the  present  analysis  from 
assumptions  of  classical  response. 


In  eq.  (2.1),  more  generally,  an  integral  over  x’is  also  present;  but  for 
local  response,  the  most  general  (linear)  Kernel  e  =  t  (3?,  t  ;3,  t ')  is  proportional 
to6(3)  (t-^'),  which  results  in  eq.  (2.1). 
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For  a  stationary  medium,  with  p,T  independent  of  t,  Eq.  (2.3)  is  a 
convolution,  so  its  Fourier  transform  with  respect  to  the  time  is 

D(*,w)  =  e(p,T;u>)E(x,u>),  (2.4) 

with  an  obvious  notation;  in  particular  e(p,t;w)  is  the  Fourier  transform  of  the 
right  side,  of  Eq.  (2.2),  but  with  p,T  independent  of  t.  When  explicit  t-depen- 
dence  from  medium  response  is  present  in  £,  Eq.  (2.4)  does  not  follow  from 
Eq.  (2.5). 

Finally,  the  dc-conductivity  a  depends  upon  the  thermodynamic  state  of  the 
fluid  also,  0=o(p,T),  and 

J(x,t)  =  o (p (x, t) ,  T(x,t) (E(x,t)  (2.5) 

is  the  quasi-static  model  analogue  of  Eq.  (2.3). 

The  electromagnetic  fields  of  the  laser  wave  in  the  fluid  will  have  a  fast 
and  a  slow  time-scale  (envelope)  part,  and  equations  for  the  latter  will  be 
derived.  Preparatory  to  attendant  labors  of  this  effort  I  list  a  few  relation¬ 
ships  issuing  from  eq.  (2.3).  The  simplest  is 

VxD  =  Ve  (p,T)xE  +  e  (p,T)VxE,  (2.6a) 

op  op 

where  the  meaning  of  the  first  term  on  the  right  side  is 

Ve  (p,T)xE=/dt'Ve(p(x,t) ,T(x,t) .t-t')xE(t') .  (2.6b) 

op 

If  E  =  0  at  differentiation  of  eq.  (2.3a)  gives 


3D 

3t 


(P.T) 


at 


•E  + 


3E 

at’ 


(2.7a) 


where 


3£0P(P,T)  E  3e(p  Cx,t).TCx,t);t-t»)E(ti;) 


(2.7b) 


in  which  t-t'  is  held  fixed  for  the  derivative  of  e.  Various  further  relations 
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follow  similarly,  such  as. 


_3 

3t 


(Ve  xE) 
op 


3Ve 


at 


2Exe 

f 


+ 


Ve  x 


op 


3E 

at' 


(2.8) 


In  what  follows,  I  will  suppress  the  p,T  dependence  of  e0p  in  the  notation,  and 
I  will  drop  the  "op"  suffix  as  well. 

Finally,  I  also  need  the  inverse  of  the  integral  operator  e,  which  I  define 

by 

/\ 

E(t)  5  (e-1  •  D)(t)  =/dt’e_1(p,T;t-t')D(t').  (2.9) 


Substituting  Eq.  (2.3a),  and  denoting  t ' -dependence  by  primes  on  p  and  T, 


/\ 

E(t)  =/ dt'/dt"c‘1(p,T,t-t')e(p',T';f-t")E(t") 
whence,  since  E  is  arbitrary. 


/\ 

-1 


/dt'  e'1(p,T;t-t,)e(p,,T,;t,-t")  =  6(t-t"). 


(2.10) 


(2.11a) 


Similarly,  starting  from  D  =  eE,  one  finds 

/dt'e(p,T;t-t,)e“1(p' ,T';t'-t")  =  6(t-t").  (2.11b) 


Notice  that 

/\ 

$ dt ' (p ,T; t-t ' )e (p ,T;t ' -t")  t  6(t-t")  (2.12) 

unless  ap/at  =  0  and  3T/at  =  0;  this  parallels  the  remark  made  earlier  concerning 
Eq.  (2.4).  Notice  also  that  the  left  sides  of  Eqs.  (2.11)  are  not  convolutions. 
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Note,  finally,  that  Eq.  (2.11),  which  I  may  writ 
-1 

e  e  =  1  , 

op 

implies 

Ve  ^  .  e=  -e  .  Ve , 

since  5(t-t")  has  no  x-dependence.  In  the  same  way. 
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(3.5) 
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reads 


,2„  .  S-il 
2  2 
c  at 


4ir  _  ,  1  3Ve  _  1_  3E  A  1  3e  3B 

-  —  Vx(oE)  -  -  JpxE  - 

c 


(3.6) 


For  the  third  term  on  the  right,  I  use  Eq.  (3.2),  with  Eq.  (2.3),  to  write 


1  _  3E  _  -1/e 

-  ?£X—  =  VtXC  - 


*(i  nK-  = 

-  — *(VxB  -  irE  -  elf E); 


(3.7) 


while  for  the  first,  I  use  Eq.  (3.1),  and  have 


4tt_  t  ci  4ir_  _  4iro  3B 
Vx(aE)  =  — VoxE  -  — 

c 


(3.8) 


Substituting  these  into  Eq.  (3.6),  one  arrives  at 

(y2 "  7  ^)B  +  (V£-£’1)xVxB  -  t(«  + 


iv^a  ♦  |0xE  +  ivee1^ 


47TO  +  ffjxE 


(3.9) 


A  somewhat  less  involved  treatment  of  Eq.  (3.3),  but  along  the  same  lines,  gives 


V2  -  £-£_ 
2  2 
<  c  at, 


\  -  Kit  * 


V(V-E) 


Lqs .  (5.9)  and  (5.10)  are  exact  consequences  of  Maxwell's  equation  for  a 
nonstationary  inhomoheneous  medium  characterized  by  constitutive  relations 
defined  bvEqs.  (2.3)  and  (2.5).  The  special  case  a-*0  and  3c/3t->-0  gives  the 
field  equations. 


(3.11) 


(3.12) 


which  agrees  with 


standard 


textbook  results  for  stationary  media. ^ 
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CHAPTER  4 


SLOW  TIME-SCALE  FIELD  EQUATIONS 


I  denote  the  laser  frequency  now  by  ui,  and  write  the  fields  as 


E  =  Re(E  exp-iwt) 


(4.1) 


B  =  Re(B  exp-iu>t), 


(4.2) 


where  E  and  8  are  spatially  and  (slowly)  time-varying  complex  envelope  fields. 
In  Eqs .  (3.9)  and  (3.10)  it  is  necessary  to  evaluate  several  integral  operator 
expressions,  such  as 


~E  =  Re/dt'  3  e e"lfa)t'  E(t'), 

3t  J  at 


(4.3) 


where  I  have  suppressed  the  x-dependence  of  E.  The  right  side  of  Eq.  (4.3)  has 
the  general  form 


[e-^/d 


dxh(p,T;x)e1<1)X  E(t 


•4 


(4.4) 


where  I  have  changed  variables  with  t'=t-x.  Expanding  E  about  x=0  gives  the 
quantity  in  the  square  bracket  formally,  with  E (n) (t) =anE/3tn,  as 
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-iut 


n=o 


h(p,T;t) 


CD 

•e-iutE  ST  E^Ctl^/dte^tp.T;,) 
n=o  dw  ^ 


[• 


-iut|c,. ,,  ,  pp  n  i3E  dh(p,T;a)) 

=e  |E(t)h(p,T;u)  ♦  + 


•••]■ 


(4.5) 


where  h  is  the  Fourier  transform  of  h.  Thus,  for  Eq.  (4.3), 


h  =  h  (p  ,T;u)  S 

3  t 


(4.6) 


In  Eq.  (4.5)  the  ratio  of  the  second  term  to  the  first  is  of  order  (wt)  ,  where 
t  is  a  time-scale  for  significant  variation  of  E;  since  this  ratio  is  extremely 
small  for  the  adiabatic  conditions  of  the  problem,  only  the  first  term  has  to  be 
retained.  In  this  way,  Eq.  (4.3)  becomes,  to  very  good  approximation. 


3  2e 
3t2 


E=Re  e 


■  iwt  3  e(o)) 

3  t2 


E  , 


(4.7) 


where  I  have  written  e(w)  for  etp.Tjw)  to  simplify  notation. 
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Similarly,  the  term  containing  3  E/0  t  in  Eq.  C 3 .  10)  involves  an  expression 
of  the  form 


y*dfh(p,T;t-f)^T(e'iu)t,E(f)^ 


^ydth  (p ,  T ;  t)  elajT^- iu+|-^E 

,-i“T^(i4(n)^)“ 

n=o  \  /  dw 


hrE(t-x) 


2  2 

where  now  h  =  3  e/9  t,  while  the  3  E/9  t  -term  involves 

e"1"*  ]T  ^  /-o3-2iu£—  4Vn)  <«  •  (4. 

n=o  n’  V  3  t  /  dw 

in  which  h  =  e.  The  lowest  terms  in  Eq.  (4.8),  out  to  the  same  order  as  the 
3 2e/3  t2  -term  in  Eq.  (4.7),  include  the  n  =  1  coefficient  of  -iw,  so  that 

3  e  3  E  -  iwt  |  .  3  e  (w)  /,  .  d  \  3  e  (w)  3  E  (  ,  „ 

3trt’{Ree  [-”«  ‘(‘-“srjTt  3tJ|>  (4- 


while,  again  to  the  same  order,  Eq.  (4.9)  gives 


3  2E  I  -iwt  2  , 

e — -j  -  R  el  e  -we  (w)  fc 

9  t  / 


-  2iu(i  + 1  w1)c(w)34 


/  d  1  2d2  \  .3  2l 

v  r^r\7 
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I  neglect  terms  involving  the  second  time  derivatives  of  adiabatically 
varying  quantities  from  here  on.  Collecting  results,  Eq.  (3.10)  finally 
reduces  to 


(4.12) 


where 


e^(to)  =  e (w)  +  ^  (4tt a  +  2  |-|) 


4tt  9  q 
2  9  t' 


(4.13) 


Three  of  the  terms  in  Eq.  (3.9)  for  the  magnetic  field  involve  compositions 
of  integral  operators;  in  addition,  they  all  involve  e“*.  Each  of  these  features 
results  in  further  complications  requiring  special  attention.  The  term  containing 
VxB ,  on  the  left  side,  involves 


(Ve-e  X)  x  curl 


(t,t~) 


x  curl  e"1Uit  8  (t"). 


(4.14) 


*Eq.  (3.10)  asserts  the  vanishing  of  an  expression  of  the  form  R(e  F)  , 
where  F  is  the  left  side  of  Eq.  (4.12).  Since  F  is  sensibly  constant  for 
variations  of  t  on  the  scale  of  u>_1,  this  condition  evaluated  at  t  =  0  and 
t  =  tt/2uj  implies  the  separate  vanishing  of  real  and  imaginary  parts  of  F, 
whence  follows  Eq.  (4.12). 
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where 


Ve  e_1(t,t")  =  /  dt'Ve  (p,T;t-t')  e_1(p',T';t'-t") 


Proceeding  as  before,  with  the  variable  change  t~=t-x,  and  with  the  Taylor 
series  expansion  for  8(t-x), 


l_Ve  e  *)x  curl  B  =  Re  e 


^dT(Vc  e1)  (t. 


t-x  )e  x  curl  8(t-x) 


- iwt  S  '  in  dn  /  f,  iwx  -1 

=  Re  e  (7dTe  (Ve  e 

n=0  n<  dJV 


x  curl  ftl 


■ 
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Substituting  Eq.  (4.15),  with  some  arranging  of  the  arguments  of  the  exponential, 
and  of  p ",  T", 

Jdx  eiaJT(Ve-e_1)  (t.t-r)  =  /dtjdt"  eiu)(i>°  7e(p,T;t-t") 


?iw(t"-t  +t)e-l  ,T(t-(t-t')) ;  t'-t+x) 

/>  iait.  r  iwt 

=  /dtj  e  iVe  (p.Tjt^d^  e 

/s  ‘  V 

e"1/p(t-t1),T(t-t1);T1J 


(4.17) 


/s 

where  in  the  last  line  I  have  performed  a  Taylor  expansion  again,  this  time  of  e_1 
about  t^  =  0.  If  I  now  adopt  the  notation 


iux 


1 


e  ^oi). 


(4.18) 
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t 


N 


Eq.  (4.17)  takes  the  more  economical  form 


J'i dxe^1  (Ve.e  ^)  (t,t-r) 


.  n  ,n„  f  .  n  -1 ,  . 
l  d  7e  (u)  d  e  (a)) 


n-'  j  n 

n=0  dw 


(4.19) 


3  t 


As  noted  earlier,  Eqs.  (2.11)  do  not  have  the  form  of  convolutions,  a 
consequence  of  which  is  that  e“l(u))  j-  l/e(u)),  as  would  be  the  case  for  stationary 
medium  response.  Starting  from  Eq .  (2.11a),  with  t"  =  t-T,  the  analysis  leading 
to  Eq.  (4.19)  can  be  applied,  mutatis  mutandis,  to  give 


•  n  ,n  — 1.  .  .  n  .  . 

d  e  (cj)  3  e  (oi) 


(4.20) 


dco 


3  t 


from  which  a  formal  expression  of  e-l(u)  can  be  constructed  recursively.  In  the 
lowest  order,  adiabatic  approximation,  Eq.  (4.20)  reads 


e  *(aj)e(w)  =  1  -  i 


■  dg  ^  ((d)  3  e  (iu) 


(4.21) 


du) 


*  3  t 


so  that 


e  1  (w) 


1_ _ 

e(w) 


-  l 


e  (w)  dw 


(  1  \  3  e(cj) 

\e(m)/H 


(4.22) 


which  is  sufficient  for  purposes  of  the  present  problem.  Eq.  (4.20)  and 
Eq.  (4.22)  both  reduce  to  the  stationary  medium  result  in  the  limit  of  e(w) 
independent  of  t. 


21 


NSWC  TR  83-100 


combining  Eq.  (4.19)  and  (4.22),  one  has 


fd tel“T  ‘  j^y  -  ]  ,eM 

♦  13T(rW)sr,EM'  ( 


(4.23) 


for  use  in  Eq.  (4.16);  substituting,  one  arrives,  after  a  little  algebra,  at 


(VEE'1) 


x  curl  B  ±  Re\e 


-icot  1 


Ve(w)x  curl  8(t) 


-i  ^.og  g.Cco^  d_  v  log  x  curl  g  (t) 


+  i  ^Yv  log  e(w)j  x  curl 


(4.24) 


That  leaves  the  last  two  of  the  four  terms  on  the  right  side  of  Eq.  (3.9) 
The  first  of  these  is  similar  to  (4.16), 


^Ve-c  ^  xaE  =  Rej^e  J" dx  Ve  •  e  *  (t,t-x)e^aJT  xa(t-x)E(t-x)^j 

=  Re[e’la)ty  rr  /*dxelu,T(Ve  .e-1)  (t,t-x)\ 
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whence,  with  the  aid  of  Eqs.  (4.19)  and  (4.22), 


^Ve-e^xaE  =  Re|e  ^-jjVe  (uj)xE 

-ioLlo|F£Md_(71°g  e(a)^xE  +  il-^vlog  c(^))x|T(aE)j|  (4. 


26) 


The  remaining  term  contains  a  product  of  three  operators;  but  it  is  also  already 
of  first  order  smallness  in  the  adiabatic  regime.  One  has 


Ve  ■  e 


■1  3e  „  „  f  -iajt  /*,  iur/^  -1  3q  . 

•  — xE  =  Re  e  J  die  IVe-e  •  —  l(t,t-T) 


xE(t) 


] 


(4.27) 


with 


3  e(u>) 

3  t 


1 

e(w) 


Ve  (u) 


3  e  (co) 
3  t 


(4.28) 
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v. 


vlvlv^vlvlv 


»  *  * 
rjl  \Jk 
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so  that 


Ve  •  e  *  •  x  E  =  Re 
d  t 


r  -iwt  3  e ( 
L6  3  t 


Vlog  e(o>)  xE  (t) 


(4.29) 


The  equation  to  be  satisfied  by  B(t),  analogous  to  Eq.  (4.12)  for  E(t),  c 
be  written  down  directly  after  substituting  Eq.  (4.2)  into  Eq.  (3.9)  by  making 
use  of  the  foregoing  results,  viz.  Eqs.  (4.11),  (4.24),  (4.10),  Eq.  (4.5)  with 
h=  3e(w)/3t,  and  Eqs.  (4.26)  and  (4.29).  That  equation  is 


where 


[v2  +  ^  e CB)  (oo)j  B  +  ^1  -  i  ^Vlog  e (to)  x  curl  B 

.  d  /  .  ..  ,\  ,  38  2ito/.  u  d  \  r  ,  3B 

+  1  d^V1°g  e(^jX  curl  3~  +  ~2  y-  +  2  d^)E(u,)'  3~T 


=  a(e,o)  x  E  +  b(e,o)  x 


(4.30) 


e  ^  (w)  =  e  (w)  +  —  (  4ttct  +  | 
0)  \  d 


e  (o>)  ^ 


(4.31a) 


a(e,o)  = 


4tt„  .  a  ( 
— Vo+  —  ( 
c  c  \ 


,  .  3  log  e(w)  d  \„,  ,  •> 

1  ♦  i  —ft - J^jviog  e(«) 


.  3  o  d  nl  ,  A  £  9  _  ,  f 

+  1  3~  dwV1°8  £  (u)  '  c3~V  108  E(W) 


(4.31b) 


b  (e  ,a)  =  i  ~  V  log  e(w) . 


(4.31c) 


Eqs.  (4.12)  and  (4.30)  are  the  slow  time-scale  equations  of  the  wave.  No 


approximations  have  been  made  regarding  spatially  varying  properties  of  the  medium. 


NSWC  TR  83-100 


I 


»  .  •  .  *t.v> 


CHAPTER  5 


ONE-DIMENSIONAL  MODEL 


1  assume  a  one-dimensional ly  stratified  fluid  now,  for  which  e  depends 
significantly  upon  only  one  of  the  three  spatial  variables,  z,  so  that 


„  “  3  e 

VE  =  2  3T  ’ 


(S.l) 


where  z  is  a  unit-vector  along  the  positive  z  -  axis  (downward  into  the  fluid,  let 

us  say.)  Eqs.  (4.12)  and  (4.30)  are  consistent  in  this  case  with  vanishing 
z  -  components  of  the  field  amplitudes.  Assuming  Ez=  0,  one  can  derive  V-c  =  0 
from  V-V  =  0;  one  has 


V-  (eg)  =  Ve-E  +  e  V-g  =  0 


(5.2) 


so  that 


V-E  =  -e  Ave-E  =0,  E  =  0. 

z 


From  Eq.  (5.3)  and  Eq.  (4.1),  Eq.  (4.12)  then  becomes 


(5.3) 


P  •  i 


e(E)(aj)  Ex  +  - 


1  + 


\ 

f  ^  JT  - 


(5.4) 


with  a  similar  equation  for  Ey  .  I  consider  here  only  the  case  that  the  direction  of 
E  is  constant  and  uniform,  hcwever,  which  I  take  now  to  be  the  x-axis, 


E  =  E  x,  E  =  0, 
x  y 


(5.5) 


wt 
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ThenEq.  (4.30)  supports  solutions, 


8  =  By  y  ,  y  =  zxi. 


(5. 


From  8=0, 
z  1 


V  log  e  (w)  x  curl  8  = 


9  log  e(oj)  y 
3Z  92  ’ 


(5. 


with  a  similar  equation  for  the  curl  (3  8/3  t)  -  term  in  (4.30).  The  equation  that 
results  for  8  is 

y 


.  9  log  e(u)  d  \  3  log  e  (u)  3 _  uj_ 

1  3  t  da)  J  3  z  3  z  +  2 


By 


-  d_  / 3 _ log  c(oj)\  3  2 By  2  ioi 

du  \  3  z  /  3  z3  t 


3By 

3t 


=  a(e,a)  +  b  (e,a)  E  . 


(5. 


Eq.  (5.4)  for  Ex  is  readily  interpreted.  Its  time- independent  form,  with 
also  a  =  0,  is 


(5. 


by  Eq.  (4.13).  This  is  the  usual  time- independent  "Schrodinger"  equation  for 
propagation  in  an  inhomogeneous  medium,  with  e  playing  the  role  of  (optical) 
potential.  The  eikonal,  or  WKB  approximation  to  Eq.  (5.9),  for  the  case 
V2  ->  32/3z2,  is 


E 


'xfl 


exp 


( -K 


z 

dz 


«/e(p(z'),T(z');u>) 


) 


(5.10) 
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where  k  £  u/c  and  Exo  is  a  constant,  and  where  the  positive  square  root  (forward 
propagating,  damped  solution)  is  understood.  When  the  WKB  approximation  is  invalid, 
which  is  the  case  of  significant  z-dependence  of  e  on  the  scale  of  c/w  =  2  ■£, 

the  solution,  Eq.  (5.10),  no  longer  holds,  and  Eq.  (5.9)  has  to  be  solved  explicitly. 
Finally,  in  Eq.  (5.4),  the  time  derivative  term  reflects  a  wave  dynamical  effect 
associated  with  medium  non-stationarity . 

Eq.  (5.8)  for  By  is  more  complicated,  a  situation  also  true  in  the  time- 
independent,  stationary  medium  limit.  In  the  latter  case,  with  also  o  =  0  again, 

Eq.  (5.8)  becomes 


( 


-> 

V" 


3  log  £  (cn) 
3  z 


(a.)  3 


0. 


(5.11! 


The  substitution. 


6  =  Je(co)  6  ,  (5.12) 

y  1  y 


eliminates  the  3 /3z  terms  from  (5.11)  giving 


('*  'i  ™)*y  ■  O' 


(5.13) 


the  same  form  as  Eq.  (5.9),  but  with  now 


G  («*)) 


£  (w) 


ci  r  1  3_2_i 

oj2  L2  3  z 


log  e(w) 
2 


—  (  9  log  £  (wi) 
4  \  3  z 


(5.14) 


This  recovers  the  text-book  result  again, ^  including  the  equality  e*  -  e  in  the 
geometrical  optics  approximation.  For  nonstationary  media,  the  substitution  of 


Landau,  L.  D. ,  and  Lifshitz,  E.  M. ,  Electrodynamics  of  Continuous  Media, 
(Reading,  Massachusetts:  Addison-Wesley ,  1960,  p.  284ff. 
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(.5.12)  does  not  accomplish  the  same  effect  owing  to  the  presence  of  the  cross¬ 
derivative,  3  2/3  t3 z -term  in  eq.  (5.8).  Evidently  also,  this  term  reflects  an 
additional  coupling  effect  between  medium  inhomogeneity  and  nonstationarity. 

Eqs.  (5.4)  and  (5.8)  have  a  somewhat  different  charactor  due  to  the 
cross-derivative  term  in  the  latter.  But  the  former  is  closely  related  to  the 
equation  for  propagation  in  a  medium  with  spatially  nonuniform  gain,  viz. 


(v?  -  2  ik  *-).  .  -i 


■  ikGe , 


(5.15) 


2  2  2  2  2 

where  V  =  3  /Ox  +  3  /3y  ,  e  =  e(x,y,z)  is  the  envelope  emplitude  and  G  =  G(x,y,z) 
is  the  gain  function.  Eq.  (5.15)  has  recently  been  solved  analytically,  for 
specified  G,  by  B.  N.  Perry,  et  al.,17  and  techniques  introduced  by  these  workers 
should  be  useful  for  limiting  model  forms  of  eq.  (5.4). 

Eqs.  (5.4)  and  (5.8)  both  should  be  amenable  to  treatment  in  a  WKB  approxi¬ 
mation  at  least,  which  will  shed  further  light  on  the  physical  nature  of  essential 
effects  associated  with  medium  nonstationary.  I  plan  to  do  this  later. 


*The  substitution  that  does  it  is  g  =  ufc>  ,  with  u  being  a  solution  of 
,3u  .  Su  n  y  y 

2  3—  *  bl  ST  aiu  '  °- 

where  a^  is  the  coefficient  of  3  8y/)z  and  b^  the  coefficient  of  3  By/)  z3t. 


17 


Perry,  B.  N.f  Rabinowitz,  P.,  and  Newstein,  M. ,  Phys.  Rev.  Lett.,  Vol.  49,  1982, 
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CHAPTER  6 


SUMMARY  AND  CONCLUSION 


I  have  derived  slow  time-scale  source-free  Maxwell  equations  for  a  nonstation 
ary,  inhomogeneous  medium,  assuming  a  high  frequency  carrier  and  slowly  varying 
response  functions  for  the  medium.  The  electric  and  magnetic  complex  field 
amplitude  equations  are  Eq.  (4.12)  and  (4.30). 

In  nearly  all  work  done  on  the  physics  of  the  coupling  of  laser  energy  into 
fluid,  the  complex  dielectric  function  e(w)  is  assumed  spatially  uniform  and 
constant.  In  fact,  it  is  a  function  of  fluid  temperature  and  density,  which  in 
many  practical  cases  are  functions  of  position,  with  steep  spatial  gradients, 
and  of  time.  Where  the  assumptions  of  homogeneity  and  stationarity  have  not  been 
made,  as  in  Ref.  8,  a  WKB  approximation  has  been  used.  The  equations  derived  in 
the  present  report  are  designed  to  permit  exploration  of  effects  for  energy 
coupling  dynamics  of  these  factors. 
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